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Abstract 

We apply the Lagrangian approach to study surface transport and 
mixing in the ocean. New tools have been developed to track the 
motion of water masses, their origin and fate and to quantify trans- 
port and mixing. To illustrate the methods used we compute the 
Lagrangian synoptic maps a comparatively small marine bay, the Pe- 
ter the Great Bay in the Japan Sea near Vladivostok city (Russia), 
and in a comparatively large region in the North Pacific, the Kuroshio 
Extension system. In the first case we use velocity data from a Japan 
Sea circulation numerical model and in the second one the velocity 
data are derived from satellite altimeter measurements of anomalies 
of the sea height distributed by AVISO. 

Keywords: Mixing; Eddy; Lagrangian synoptic map; Marine bay; Kuroshio 
Extension. 



1 Introduction 



The ocean presents a variety of dynamical phenomena with different space 
scales ranging from millemeters to a few thousand of kilometers. Despite of 
that, large-scale coherent structures are easlily visible, say, at satellite images 
of the sea color and surface temperature and can be identified by means of 
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in-situ measurements. The striking examples are the major western bound- 
ary oceanic currents, the Gulf Stream in the Atlantic and the Kuroshio in the 
Pacific. They are "rivers" with the warm water in the ocean with the width 
on the order of 100-200 km and the maximal speed of current at the surface 
of 2 m/s. Such currents separate waters with different physical, chemical and 
biological characteristics. The other examples are mesoscale (with the size 
of a few hundred of kilometers) and submesoscale (a few tenth of kilometers) 
eddies that can transport water over hundreds and even thousands of kilo- 
meters and can survive for months before breaking down. Being coherent 
features, they do not contain the same waters but exchange them with the 
surrounding ocean, the process known as mixing. 

Lagrangian and dynamical systems methods have been developed to study 
large-scale transport and mixing in the ocean jU [5j [3j |2l [6J [121 E] ■ The 
main purposes of those studies are to track the fluid motion, to elucidate and 
quantify transport and mixing processes. Simply speaking, we would like to 
know where these or those waters come from, what is their fate and how 
they mix in this or that region. In the Lagrangian approach one integrates 
trajectories for a large number of synthetic particles advected by an Eulerian 
velocity field 

Tt (1) 

The velocity field, v(f,t), is supposed to be known analytically, numerically 
or estimated from satellite altimetry. While in the Eulerian approach we get 
frozen snapshots of data, Lagrangian diagnostics enable to quantify spatio- 
time variability of the velocity field. It has been established theoretically and 
experimentally that even a simple deterministic velocity field may cause prac- 
tically unpredictable particle trajectories, the phenomenon known as chaotic 
advection [HI E]. The real oceanic flows are not, of course, deterministic 
and regular, but if the Eulerian correlation time is large as compared to the 
Lagrangian one, the problem may be treated in the framework of chaotic 
advection concept. 

It is important to separate chaotic and turbulent mixing in the ocean. 
The process of chaotic advection provides transport and mixing with the 
characteristic scales on the order of a few tenths or even hundreds of kilome- 
ters, whereas turbulence works at smaller scales. At a comparatively large 
scale, turbulent mixing is homogeneous whereas the chaotic one is not. Typ- 
ical patterns of chaotic advection consist of large-scale convoluted curves 
visible in some surface-temperature and color satellite images. The effect 
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of turbulent mixing is in small-scale fluctuations superimposed on the large- 
scale convoluted curves. If the velocity field on comparatively large scales 
is quasicoherent in space and quasiregular in time but the motion of tracers 
is mainly irregular, one deals with chaotic mixing. Turbulent mixing means 
that the velocity field is irregular in space and time at the same scales at 
which the tracer's motion is irregular. 

In this paper we report on our recent results on developing Lagrangian 
tools to monitor surface transport and mixing in the ocean. We propose with 
this aim new Lagrangian criteria that enable to track and quantify the water 
exchange processes and reveal the underlying physical mechanisms. As an 
output, we compute different Lagrangian synoptic maps of the regions under 
study for a given period of year and analyze them. The methos is illustrated 
with a comparatively small marine bay, the Peter the Great Bay in the Japan 
Sea near Vladivostok (Russia), and a comparatively large region in the North 
Pacific, the Kuroshio Extension system. In the first case we have used velocity 
data from a Japan Sea eddy-resolved circulation numerical model with the 
fine resolution of 2.5 km, in the second one — satellite altimetric velocity 
data with the coarse resolution of the order of 35 km. 



2 Lagrangian and dynamical systems meth- 
ods to study transport and mixing in the 
ocean 

Motion of a fluid particle in a two-dimensional flow is the trajectory of a 
dynamical system with given initial conditions governed by the velocity field 
computed either by solving the corresponding master equations or as the 
output of a numerical ocean model or derived from a measurement 

^ = u(x, y,t), ^ = v(x, y, t), (2) 

where (x, y) is the location of the particle, u and v are the zonal and merid- 
ional components of its velocity. Even if the Eulerian velocity field is fully 
deterministic, the particle's trajectories may be very complicated and prac- 
tically unpredictable. It means that a distance between two initially nearby 
particles grows exponentially in time 

||<5r(0IIHIM0)||e A *, (3) 
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Tp — y^ 11 nh?^i (4) 



where A is a positive number, known as the Lyapunov exponent, which char- 
acterizes asymptotically the average rate of the particle dispersion, and || • || 
is a norm of the vector r = (x, y). It immediately follows from ([3]) that we are 
unable to forecast the fate of the particles beyond the so-called predictability 
horizon 

A ||A(0)|| 

where ||A|| is the confidence interval of the particle location and ||A(0)|| is a 
practically inevitable inaccuracy in specifying the initial location. The deter- 
ministic dynamical system ([2]) with a positive maximal Lyapunov exponent 
for almost all vectors 5r(0) (in the sense of nonzero measure) is called chaotic. 
It should be stressed that the dependence of the predictability horizon T p on 
the lack of our knowledge of exact location is logarithmic, i. e., it is much 
weaker than on the measure of dynamical instability quantified by A. Simply 
speaking, with any reasonable degree of accuracy on specifying initial con- 
ditions there is a time interval beyond which the forecast is impossible, and 
that time may be rather short for chaotic systems. 

Since the phase plane of the two-dimensional dynamical system (T5]) is 
the physical space for fluid particles, many abstract mathematical objects 
from dynamical systems theory (stationary points, KAM tori, stable and 
unstable manifolds, periodic and chaotic orbits, etc.) are material surfaces, 
curves and points in fluid flows. It is well known that besides "trivial" el- 
liptic fixed points, the motion around which is stable, there are hyperbolic 
fixed points which organize fluid motion in their neighbourhood in a specific 
way. In a steady flow the hyperbolic points are typically connected by the 
separatrices which are their stable and unstable invariant manifolds. In a 
time-periodic flow the hyperbolic points are replaced by the corresponding 
hyperbolic trajectories with associated invariant manifolds which in general 
intersect transversally resulting in a complex manifold structure known as a 
heteroclinic tangle. The fluid motion in these regions is so complicated that 
it may be strictly called chaotic, the phenomenon known as chaotic advection 
[HI [5]. Adjacent fluid particles in such tangles rapidly diverge providing very 
effective mechanism for mixing. 

Stable and unstable manifolds are important organizing structures in the 
flow because they attract and repel fluid particles (not belonging to them) at 
an exponential rate and partition the flow into regions with different types of 
motion. Invariant manifold in a two-dimensional flow is a material line, i. e., 
it is composed of the same fluid particles in course of time. By definition 
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stable (W s ) and unstable (W u ) manifolds of a hyperbolic trajectory j(t) are 
material lines consisting of a set of points through which at time moment t 
pass trajectories asymptotical to 7(i) at t — Y oo (W s ) and t — > — oo (WJ. 
They are complicated curves infinite in time and space that act as boundaries 
to fluid transport. 

The real oceanic flows are not, of course, strictly time-periodic. However, 
in aperiodic flows there exist under some mild conditions hyperbolic points 
and trajectories of a transient nature. In aperiodic flows it is possible to iden- 
tify aperiodically moving hyperbolic points with stable and unstable effective 
manifolds [U [3] . Unlike the manifolds in steady and periodic flows, defined 
in the infinite time limit, the "effective" manifolds of aperiodic hyperbolic 
trajectories have a finite lifetime. The point is that they play the same role 
in organizing oceanic flows as do invariant manifolds in simpler flows. The 
effective manifolds in course of their life undergo stretching and folding at 
progressively small scales and intersect each other in the homoclinic points 
in the vicinity of which fluid particles move chaotically. Trajectories of ini- 
tially close fluid particles diverge rapidly in these regions, and particles from 
other regions appear there. It is the mechanism for effective transport and 
mixing of water masses in the ocean. Moreover, stable and unstable effective 
manifolds constitute Lagrangian transport barriers between different regions 
because they are material invariant curves that cannot be crossed by purely 
advective processes. 

The stable and unstable manifolds of influencial hyperbolic trajectories 
are so important because (1) they form a kind of a sceleton in oceanic flows, 
(2) they divide a flow in dynamically different regions, (3) they are in charge 
of forming an inhomogeneous mixing with spirals, filaments and intrusions, 
(4) they are transport barriers separating water masses with different char- 
acteristics. Stable manifolds act as repellers for surrounding waters but un- 
stable ones are attractors. That is why unstable manifolds may be rich in 
nutrients being oceanic "dining rooms" . 

There is a quantity, the finite-time Lyapunov exponents (FTLE), that 
enables to detect and visualize stable and unstable manifolds in complex 
velocity fields. The FTLE is the finite-time average of the maximal separation 
rate for a pair of neighbouring advected particles which is given by [9] 

X(r(t)) = -lna(G(t)), (5) 

T 

where r is an integration time, a(G(t)) the largest singular value of the evo- 
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lution matrix for linearized advection equations. Scalar field of the FTLE is 
Eulerian but the very quantity is a Lagrangian one that measures an inte- 
grated separation between trajectories. Ridges (curves of the local maxima) 
of the FTLE field visualize stable manifolds when integrating advection equa- 
tions forward in time and unstable ones when integrating them backward in 
time. 

3 Transport and mixing in marine bays 

When studing transport and mixing in marine bays, it is important to know 
which waters enter the bay under study, which ones quit the bay, by which 
transport corridors they do that and how the different waters mix in the 
bay interior. The Lagrangian approach, allowing to compute the origin and 
fate of different waters, is the most suitable for that. Transport and mixing 
in marine bays is more inhomogeneous as compared with those processes 
in open basins because of a complicated structure of currents and eddies of 
different scales, strong tides and presence of river estuaries. In this section we 
apply Largangian tools to characrerize horizontal subsurface transport and 
mixing in the Peter the Great Bay near Vladivostok city (Russia). That is 
the largest bay in the Japan Sea with a few shallow-water smaller bays and 
estuaries of three major rivers with a wide shelf and steep continental slope. 
The water exchange between the bay and the open sea is governed mainly by 
a cyclonic circulation over the deep central basin and the Primorskoye current 
flowing to the southwest along the continental slope of the Primorsky Krai 
(Russia). We have used velocity data from the MHI ocean circulation model 
[TU] which is a set of 3D primitive equations in Z-coordinate system with 10 
quasi-isopycnal layers and the resolution of 2.5 km. 

To characterize the water exchange between the Peter the Great Bay and 
the open sea we compute the FTLE map and the exit-time map (Fig. [Q. 
A large number of synthetic particles have been uniformly distributed over 
the region with [130°12' : 133°12'] E and [41°42' : 43°19'] N. In Fig. Q}l we 
compute the FTLE, A, by the method proposed in Ref. [9]. The advection 
equations (J2J) have been integrated forward in time for 54 days in the August 
and September of a typical year. The gray shades code the magnitude of A. 
The value A = 0.085, at which the distance between neighbouring particles 
increases in 100 times, is chosed to be a threshold. The regions with A < 0.085 
are supposed to be regular, the ones with A > 0.085 — chaotic. The black 
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ridges with A 3> 0.085 visualize stable manifolds of influencial hyperbolic 
trajectories in the region. Spiral-like structures reveal eddies of different 
scales, the white and light-grey zones are the stagnation regions or shear 
currents. The sandwich-like structures are signs of the most intense mixing. 
The synoptic Lyapunov map in Fig. [T^, shows the scalar filed of this quantity 
in geographic coordinates which are initial positions of the synthetic particles. 
This map along with the Lyapunov map, computed backward in time (not 
shown), demonstrates with a high resolution the complicated character of 
transport and mixing in the Peter the Great Bay. 

The exit-time map is shown in Fig. [lb. The color in the map codes 
the time, T, particles (initially distributed over the same region) need to 
reach the open sea or the coastline. In fact, we compute the trajectories 
till they reach the 3 km band along the coastline. The white wide band 
along the coast in Fig. [lb demonstrates the Primorskoye current along which 
particles quickly leave the bay to the southwest. The large white corridor 
in the central part of the region selected separates the Peter the Great Bay 
from the open sea. Black color marks the particles that did not leave the 
bay for the computation time, 54 days. The stagnation zones are situted, as 
expected, in the smaller bays, the Amursky and the Ussyrisky ones, which 
are visible as black spots on the both sides of the peninsula in the north. 
The exit-time map reveals the complicated process of chaotic mixing in the 
central part of the bay with the spiral-like anticyclonic eddy (with the center 
at 132°45' E and 42°40' E) and gives a valuable information about origin and 
fate of waters. 

To get an information about the character of motion of different waters, 
their drift, rotation and oscillation, we compute the new Lagrangian synoptic 
maps: rotation and mixing maps, transport and visitor maps. We compute 
for a large number of particles the number of cyclonic, r) c , and anticyclonic, rj a , 
rotations and their difference i]. The typical kinds of particle's motion are the 
following: 1) simple drift or linear displacement if r] c , r] a , \r)\ < r] cr = 5, where 
r] cr is a threshold value of the rotation number; 2) rotation, if \rj\ > r] cr = 5; 
3) oscillation, if r] c , r\ a > r] cr = 5 but \rj\ < r] CT . In the rotation map in Fig. [2^ 
white and black colors mean cyclonic, r] c , and anticyclonic, r] a , rotations, 
respectively, computed for the same period of time, 54 days. Grey color 
codes the particle with predominant displacements or oscillations. The map 
demonstrates clearly the same spiral-like anticyclonic eddy as in Fig. [1] and 
the large-scale filaments with foldings typical to chaotic advection in the 
ocean. 
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Figure 1: (a) The Lyapunov map in the Peter the Great Bay and the sur- 
rounding region of the Japan Sea. (b) The exit-time map in the same region. 

To characterize the chaotic mixing more clearly we compute along with 
the rotation numbers the FTLE A. If A > A cr = 0.85 and 7] a > r] cr = 5, we will 
speak about unstable rotations in the corresponding region. If A > A cr = 0.85 
but r) a < r] CI = 5 one deals with unstable linear displacement of the corre- 
sponding particles. The mixing map in Fig. [2b shows by color regions with 
different dynamical properties specified by the rotation numbers and the 
maximal Lyapunov exponent. White color marks the regions with regular 
oscillations and/or predominant displacements. The spots of particles, placed 
in those regions, move as whole being deformed slightly. The white grey color 
- the regions with unstable displacements which are peripheries of the an- 
ticyclonic eddies and their filaments. The spots, placed in those regions, are 
elongated strongly. The dark grey color — the regions with unstable oscil- 
lator motion with the particles rotating for 54 days in the cyclonic and then 
in the anticyclonic directions. The black color corresponds to the unstable 
rotation that manifests itself in narrow filaments and spiral-like structures in 
anticyclones. 

In order to find frontal zones and transport pathways we propose to com- 
pute the transport maps showing the final positions of particles when in- 
tegrating the advection equations (J2]) forward and backward in time (see 
Figs. |3k and|3b, respectively). In other words, the equations © have been 
solved for each of the million particles initially distributed over the region 
selected for 54 days forward and backward in time. In the first case we get 
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Figure 2: (a) Rotation and (b) mixing maps in the Peter the Great Bay. 

the particle's fate map (Fig. [3^) with the black (white) particles leaving the 
bay through the eastern (western) border. The grey particles are those that 
did not leave the bay for the computation time. When integrating the equa- 
tions 02]) backward in time, we get the particle's origin map with the black 
(white) particles entering the bay through the eastern (western) border and 
the resident particles shown in grey. The frontal zone, separating the waters 
with different fate and origin, consistes of smooth, meandered and spiral-like 
fragments. 
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Figure 3: Transport maps in the Peter the Great Bay. (a) The particle's fate 
and (b) origin maps. 
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4 Transport and mixing in the Kuroshio Ex- 
tension region 

The Kuroshio Extension prolongs the Kuroshio Current when the latter sep- 
arates from the continental shelf at about 30° N. It flows eastward from this 
point as a strong unstable meandering jet constituting a front separating the 
warm subtropical and cold subpolar waters of the North Pacific Ocean. It 
is a region with one of the most intense air-sea heat exchange and the high- 
est eddy kinetic energy level strongly affecting climate. Transport of water 
masses is of cruicial importance and may cause heating and freshing of waters 
with a great impact on the weather and living organisms. 

The surface ocean currents used in this section are derived from satel- 
lite altimeter measurements of sea height (http:/ /www.aviso. oceanobs.com ). 
The velocity data covers the period from 1992 to 2011 with weekly data on 
a 1/3° Mercator grid. In our study we focus on the region between 30° and 
45° N and between 130° and 165° E. Bicubical spatial interpolation and third 
order Lagrangian polinomials in time have been used to provide accurate nu- 
merical results. Lagrangian synoptic maps, manifolds and chaotic advection 
structures in general are determined by the large-scale advection field, which 
is appropriately captured by altimetry. Thus, computation of particle's tra- 
jectories statistically is not especially sensitive to imperfections of the velocity 
field caused by the interpolation and measurement imperfections. 

In Fig. we demonstrate the displacement map for the region computed 
for 45 days after the beginning of the incident at the Fukushima Daiichi 
nuclear power plant. The shades of gray depict the magnitude of the dis- 
placement of a tracer D = {xf — x ) 2 + (jjf — yo) 2 , from its initial position, 
(x ,y ), to a final one (xf,yf). The Kuroshio Current is well pronounced 
including meanders and intrusions, its extension, and mesoscale eddies. Two 
light-colored eddy patches are of particular interest. Their centers are ap- 
proximately at the latitude of the Fukushima plant and at longitudes 153° E 
(the mushroom-like dipolar eddy) and 161° E (the circular eddy), both ed- 
dies being surrounded by dark-colored necklaces having a relatively high 
magnitude of D. This pattern exemplifies the ring birth process due to the 
meandering of the Kuroshio current and subsequent detachment of eddies 
from the main jet. 

To get a picture of an "instantaneous " state of the region we show in 
Fig. Hb the surface velocity field computed on the fixed day, 1 January, 2010. 
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Figure 4: (a) Displacement map in the Kuroshio Extension region. The color 
codes the magnitude of displacement D in minutes. The Fukushima Daiichi 
nuclear power plant is marked with the sign of radioactivity, (b) Velocity 
field of the region on 1 January, 2010. Circles and crosses are instantaneous 
elliptic and hyperbolic points, respectively. 

The main meandering jet is depicted by the black arrows corresponding to 
comparatively large velocities. The grey arrows around the jet with smaller 
velocities reveal a number of cyclonic and anticyclonic eddies on both flanks 
of the jet. We have computed instantaneous elliptic and hyperbolic points of 
the flow and showed them by circles and crosses, respectively. The elliptic 
points are situated mainly in the centers of the eddies, whereas the hyper- 
bolic ones are in the regions between the eddies with different polarity and in 
the periphery of isolated eddies. The hyperbolic points are especially impor- 
tant because they may be connected by instantaneous stable and unstable 
manifolds dividing the flow into regions with cardinally different dynamics. 

We plan to show in this section that the Lagrangian diagnostics is well 
suitable to describe the mesoscale and submesoclace features of the complex 
picture of mixing in the Kuroshio Extension region. The altimetric velocity 
data we used covers the period from 1 January to 3 July, 2010. We focus on 
a vortex pair on the jet's southern flank, consisting of the anticyclone (AC) 
and cyclone (C). The pair manifests itself on the Lagrangian maps in Fig. [5] 
computed for a large number of synthetic particles seeded over the region 
considered. All the maps visualize the northern hat-like AC with the axes 
of 150 and 100 km and the southern circular C with the diameter 150 km. 
In Fig. [5^ the color codes the meridional displacement, D y , of particles on 
the 60th day of integration. The spiral structure of the C is well developed 
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with the spiral untwisting counter-clockwise, whereas it is less pronounced 
for the AC with the spiral untwisting clockwise. The character of the water 
motion in the C and AC is also different and becomes evident after computing 
the number of particle's rotation around the vortex centers. It follows from 
Fig. |5b that water in the AC core circulates with approximately the same 
angular velocity, whereas this quantity decreases from the center of the C 
to its periphery (pay attention to the ring-like structure of the C). In order 
to visualize the stable manifolds of the hyperbolic trajectories around the 
vortex pair, we compute in Fig. |5b the FTLE, A, and displacements, D, of the 
particles. The shades of grey in this figure modulate different combinations 
and magnitudes of A and D with respect to some chosen "critical" values: 
A cr , corresponding to divergence of initially close particles over 100 km, and 
D cr = 100 km. The black convoluted curves in the figure between the eddies, 
around each of them and around the very pair delineate the corresponding 
W s manifolds. 




32° 33° 34° 35° 36° 37° 38° -2 -1 1 2 3 4 5 



Figure 5: Lagrangian maps of the AC C vortex pair in which the color codes: 
(a) the meridional displacement of synthetic particles, D y , on the 60th day 
of integration, (b) the number of their rotation around the vortex centers on 
the 15th day and (c) their Lyapunov exponents, A, and, D, displacements 
on the 45th day with the following legenda: white means A < A cr , D > D CI , 
grey — A < A cr , D < D CT and black — A > A cr , D > D CT . 

To give a detailed description of the structure of each eddy in the vortex 
pair we apply the method of particle's scattering elaborated in Ref. pp. We 
cross both the eddies by a material line and compute rotation number T], and 
the maximal FTLE on initial particle's latitude y . The scattering plot in 
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Fig. demonstrates that the waters in the C core really rotate with different 
angular velocities decreasing from the center to its periphery. Rotation in the 
AC core is much more homogeneous. Moreover, the waters in the C rotate in 
two times faster than in the AC. The scattering plot \(yo) in Fig. [6j3 demon- 
strates smooth segments in the cores of C and AC and irregular oscillations 
in their periphery. It simply means that the water in the cores moves more 
or less coherently whereas the motion in the eddy's peripheries is erratic due 
to numerous intersections of stable and unstable manifolds. Computation 
of the dependence of the time of exit of the particles T, belonging to the 
material line, on yo (not shown) confirms that waters prefer to quit the C 
more or less periodically by portions. Each portion is represented by a U-like 
segment of the T(y ) function which consists of a large number of particles 
with approximately the same time of exit and the same rotation number 77. 
In difference from the C, particles quit the AC core practically at the same 
time. In other words, the particles quit the C by portions along spiral-like 
transport pathways, whereas the periphery of the AC exchanges water with 
the surrounding but its core moves coherently as a whole for a time. 




32° 33° 34° 35° 36° v 32° 33° 34° 35° 36° 



Figure 6: The scattering plots for the vortex pair on the 30th day of integra- 
tion, (a) Number of times, rj, the particles rotate around the vortex centers 
vs initial particle's latitude position y , (b) the corresponding maximal FTLE 
vs y . 

In conclusion we demonstrate in Fig. [7] how frequently fluid particles, 
chosen in the cores of the C and AC, visit for 180 days different places in 
the Kuroshio Extension region. It is evident that the C was absorbed by 
the main jet in a short time and then its waters travelled eratically within 
the jet with a few excursions to its northern and southern flanks. It is 
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interesting that in course of time C waters have formed the new cyclonic 
eddy nearby (x = 155° E, y = 32° N). In contrast to the C, the AC waters 
have walked eratically on the southern flank of the jet in a restricted region 
within x = [140° : 150°] E, y = [28° : 35°] N. 
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Figure 7: Visitor maps for (a) the cyclone and (b) anticyclone show how 
frequently fluid particles from the corresponding eddy's cores visit for 180 
days different places in the Kuroshio Extension region. 



5 Conclusion 

The Lagrangian approach has been shown to be very useful to gain new in- 
formation on chaotic transport and mixing in the ocean. We have elaborated 
new Lagrangian diagnostic tools to visualize and quantify those processes: 
the time of exit of fluid particles off a selected box, their displacements, 
the number of their cyclonic and anticyclonic rotations and the number of 
times they visit different places in the region. Along with the Lyapunov 
maps, the corresponding high-resolution Lagrangian synoptic maps of those 
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quantities, computed by solving advection equations forward and backward 
in time for different periods of the year, are new diagnostic and prognostic 
products characterizing the state of the ocean. The technique developed can 
be applied to the global ocean and its basins. 

In this paper we have focused on a comparatively small marine bay, the 
Peter the Great Bay in the Japan Sea near Vladivostok (Russia), and on 
a comparatively large region in the North Pacific, the Kuroshio Extension 
system. In the bay study in summer and autumn periods, we have used the 
velocity data from a Japan Sea eddy-resolved circulation numerical model 
with the resolution of 2.5 km. It has been shown that the Lyapunov and 
exit-time maps, the rotation, mixing and transport maps allowed to quan- 
tify and specify movement of water masses, their mixing and the degree of 
its chaoticity in the bay. Those high-resolution maps allowed to visualize 
transport pathways by which waters exit and enter the bay. 

As to the Kuroshio Extension, we have used the velocity data derived from 
satellite altimeter measurements of sea height with the corresponding inter- 
polation. The main attention has been paid to study structure, transport and 
mixing of a vortex pair with strongly interacting cyclonic and anticyclonic 
eddies. Such dipoles occur frequently in that region. We have computed 
Lagrangian synoptic maps for the time of exit of particles, the number of 
changes of the sign of zonal and meridional velocities, and for other quan- 
tities. Along with the Lyapunov map, they have been shown to be able to 
reveal the vortex structure and its evolution, meso- and submesoscale fila- 
ments, repelling material lines, hyperbolic and non-hyperbolic regions in the 
sea. In particular, we have found that the eddies have a prominent spiral-like 
structure resembling the spiral patterns at satellite images in that region. 

The work was supported by the Program "Fundamental Problems of Non- 
linear Dynamics" of the Russian Academy of Sciences, by the Russian Foun- 
dation for Basic Research (projects nos. 09-05-98520 and 11-01-12057) and 
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